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Signatures of pairing and spin-orbit coupling in correlation functions of Fermi gases 
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We derive expressions for spin and density correlation functions in the (greatly enhanced) pseudo¬ 
gap phase of spin-orbit coupled Fermi superfiuids. Density-density correlation functions are found 
to be relatively insensitive to the presence of these Rashba effects. To arrive at spin-spin correlation 
functions we derive new /-sum rules, valid even in the absence of a spin conservation law. Our 
spin-spin correlation functions are shown to be fully consistent with these /-sum rules. Importantly, 
they provide a clear signature of the Rashba band-structure and separately help to establish the 
presence of a pseudogap. 


Introduction.— Spin-orbit coupling (SOC) in supercon¬ 
ductors and superfiuids is a topic of much current inter¬ 
est m- This is in large part because there is some 
hope that (particularly in the presence of a magnetic 
field) they may relate to the much sought after spinless 
Px + ipy superfluid [4]. Two communities have united 
around these issues: those working on cold Fermi super- 
fluids with intrinsic Rashba SOC mn] and those studying 
superconductivity that is proximity induced in a spin- 
orbit coupled material [7]. To achieve this ultimate goal 
it is important to establish that a given candidate for the 
Px + iPy superfluid simultaneously exhibits signatures of 
both pairing and spin orbit coupling. This would pro¬ 
vide minimal evidence for a properly engineered ultracold 
atomic gas. One therefore needs experimental signatures 
of these simultaneous effects and this provides a central 
goal for the present paper. 

Here we address the signatures of this anomalous spin- 
orbit coupled superfluid as reflected in spin-spin and 
density-density correlation functions. Our work builds 
on the observation [HHSl that in the presence of Rashba 
SOC, pairing (in the form of pseudogap effects [Til ITS] ) is 
significantly enhanced. For this reason (and because the 
correlation functions are free of the complications of col¬ 
lective mode effects) we focus here on the normal phase. 
We show how even without condensation, the frequency 
dependent spin response exhibits features which relate to 
the Rashba ring band-structure, as well as to the pres¬ 
ence of a pairing gap. In contrast, the density response is 
relatively unaffected by SOC. Previous work has focused 
on identifying SOC without [Min] or with pairing |18j 
via the one body spectral function. As compared with 
Ref. [IH], we find less subtle features in the two particle 
response. At the very least the spin correlation functions 
provide complementary and accessible (via neutrons or 
two photon Bragg scattering [T5]) information. 

Validating any theory of correlation functions requires 
satisfying important constraints [5D]. Indeed, the ab¬ 
sence of conservation laws complicates all spin transport 
in spin-orbit coupled materials. Thus, it is extremely 
important to find underlying principles for establishing 
self consistency. To address this issue, here we use the 
Heisenberg equations of motion to derive /-sum rules 
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FIG. 1. Phase diagrams for a degenerate fermi gas without 
(a) and with (b) SOC. Plotted is T* as a function of inverse 
scattering length, llkpa, indicating where pairing first sets 
in, while Tc marks the onset of condensation to the superfluid 
phase. The plot shows the weakly interacting regime with 
l/kpa < 0. As spin-orbit coupling is turned on, pairing is 
enhanced leading to a larger pseudogap region. The inset in 
(a) indicates the temperature dependences of the component 
gap parameters (defined in the text) at unitarity, l/kpa = 0. 


for the spin-spin correlation functions, which also pro¬ 
vides important constraints on our numerical calcula¬ 
tions. While a magnetic field is necessary for arriving 
at topological order, we begin by ignoring this additional 
complication. There is a substantial literature investigat¬ 
ing correlation functions in the superfluid phase (without 
pseudogap effects) which we note here . 

In this Rapid Communication we present two main re¬ 
sults. The first is a consistent derivation of spin-spin 
correlation functions in spin-orbit coupled Fermi gases. 
The second establishes qualitative experimental signa¬ 
tures reflecting separately the presence of a pairing gap 
and of SOC. Readers interested in the experimental sig¬ 
natures need only a cursory exploration of the mathe¬ 
matical derivation that precedes it. 

Background Theory.— We consider a gas of fermions 
whose single particle Hamiltonian is iL°(k) = k'^/2m — 
/i -b Act • ]i±/m for a particle of mass m, spin-orbit cou¬ 
pling momentum A, momentum k = {k^, ky,kz), in-plane 









2 


momentum kj_ = {kx,ky,Q), and vector of Pauli ma¬ 
trices a = (o-x, cTy, az). Throughout this paper we set 
h = ks = ^- To describe our spin-orbit coupled Fermi 
gas with pairing, we use a 4 x 4 inverse Nambu Green’s 
function 


g-HK) 


(G^\K) ^ A \ 
A G^\K))^ 


( 1 ) 


that acts on the spinor 4'^ = ^Ck^, Ck^, —j 

for a fermion annihilation (creation) operator Cks(c|js) 
of spin s =t, i and momentum k. In the Green’s func¬ 
tion A is a pairing gap, the 4-vector K = (ia;,k) with 
Matsubara frequency ioj, and the non-interacting inverse 
particle Green’s function is Gq^{K) = iuj — _ff°(k) and 
hole Green’s function Gq^{K) = iay[GQ^{—K)Yi(Ty — 
iu) + Ff°(k). 

This superfiuid can be studied at the mean field 
level [HHS], where one has the usual gap equation: 



Y,^t\g{K)G^{K) 


K 


( 2 ) 


where ^ Sitj Sk ^ momentum and 

Matsubara frequencies at temperature T, and where the 
many-body Green’s function G{K) is found from the in¬ 
verse of Eq. Q 

G-\K) - Go\K) = -S(iF) = -A^GoiK). (3) 

Except for the matrix structure in the above equations, 
these are the usual definitions of the fermionic Green’s 
functions and self energy associated with BCS theory; the 
gap equation also appears in the literature [8lH3j. As in 
the cold atoms literature, we will regularize the gap equa¬ 
tion by replacing the interaction strength g with the scat¬ 
tering length a through 1/g = mjAna — 

We now want to include pair fluctuations, or pseudo¬ 
gap (pg) effects, in a fashion fully consistent with both 
the ground state and the mean field equations that have 
been extensively studied in previous work [HVH3j. The 
approach we outline below was introduced in the context 
of high temperature superconductors [I1[5S], but it has 
also been applied to spin-orbit coupled superfluids [5S]. 

Pair fluctuations lower the phase transition tempera¬ 
ture Tf. relative to its mean field value denoted by T*. 
The latter is the temperature at which the pairing gap, 
determined by Eq. ([^, first becomes non zero. (Note that 
T* does not reflect a broken symmetry state and is not a 
true phase transition). As a result, a central component 
of the present theory is that the usual Thouless condi¬ 
tion (on the t-matrix) for the instability of the normal 
phase m must be modified to include a well developed 
excitation gap. 

Imposing consistency with mean field theory for the 
pairing gap obtained from Eq. ([^ leads to a modified 


Thouless condition: 


where 


^Pg(Q) - J ^ ^ oo, as g ^ 0, 


X(g) = [g{K)G^{K - Q) 

K 


( 4 ) 

( 5 ) 


In this way, above T), the square of the excitation gap 
of the non-condensed pairs, called Apg, is to be associ¬ 
ated with A^, given by Eq. Q [^. Establishing Tc, 
however, requires that we find a constraint on A^g appli¬ 
cable at and below T^- Indeed, once T < Tc, there must 
be another contribution to the self energy E(Ar) involving 
g = 0 or superconducting (sc) pairs. Then the self energy 
of Eq. ^ can be written in terms of a t-matrix (as in the 
Thouless condition) E(Ar) = — t(g)Go(Ar — g), pS] : 

A2 

^pg(g) +4c(g), where 4c(g) = — ^S{Q) represents the 
condensate. As a consequence E(Ar) = Esc(A')-|-Epg(Ar). 
Because of Eq. Q, the quantity tpg(g) is strongly peaked 
at g = 0 and Epg(Ar) « A^^Gq{K). 

This relation implies that A^g = — 
equivalently, the pseudogap is approximated as a ther¬ 
mal gas of composite bosons. The transition tempera¬ 
ture Tc is determined as the temperature at which this 
value of Apg, when the transition is approached from be¬ 
low, intersects with the mean field gap equation, which 
corresponds to Apg above Tc in the normal state. In 
this way, in^the ordered phase the total self energy 
E(Ar) = A^Go(Ar), where A^ = A^^, -I- Apg. In the in¬ 
set of Eig. [^a), we plot the temperature dependence for 
both Asc and Apg, as well as the total excitation gap A. 

We emphasize that a central distinguishing feature of 
the present approach is that Tc is determined in the 
presence of a well developed gap at Tc. This contrasts 
with the scheme of Nozieres and Schmitt-Rink [29] and 
is similarly different from path integral-collective mode 
schemes [301 . The latter introduce Goldstone bosons, 
but importantly these do not renormalize the mean field 
transition temperature which remains at T *. 

Figure shows the calculated phase diagram, plotting 
Tc and T*, in the absence (a) and presence (b) of Rashba 
SOG. The latter is in reasonable agreement with the re¬ 
sults of Ref. [35]. We restrict these plots to the weak 
pairing side of resonance. Here we observe a greatly en¬ 
hanced pseudogap regime denoted by an enhancement of 
T* without significant enhancement of Tc- The behav¬ 
ior of T* has been attributed to the enhancement of the 
pairing attraction [HHia, due to an increased density of 
states near the minimum of the Rashba ring. Since Tc is 
obtained in the presence of a gap at Tc, stronger pairing 
(reflected in T*) is offset by an increasingly gapped den¬ 
sity of states. This leads to a relatively constant Tc as a 
function of interaction strength. 

Density/current and f-sum rules.— To characterize the 
anomalous, normal, and superfluid phases in more detail, 
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FIG. 2. The spin-spin response function as a function of frequency. Figure [^a) corresponds to a pseudogap phase without 
SOC. The peak at lower lu reflects thermally excited fermions while the second peak is at to comparable to the pseudogap 
energy scale, where pairs are now broken. Since A = 0, the spin-spin and density-density correlation functions are the same. 
In Fig. I^b) the pseudogap is set to zero with fixed A = 1.2kF- The lower energy intra-helicity band peak is enhanced and one 
sees an SOC-related peak. As shown in the inset, the low-energy threshold reflects the onset of inter-helicity band transitions, 
while the high-energy endpoint occurs when these transitions are no longer possible. In Fig. j^c), in the presence of both a 
pseudogap and SOC, one sees a combination of the effects in the previous two panels. Spectral weight is transferred to yield 
a larger SOC peak, reflecting the gapping of the low energy contributions. All quantities are measured relative to the Fermi 
energy, Ef, or Fermi momentum fe_F. 


we investigate both the density/current and spin corre¬ 
lation functions, considering the former first. For sys¬ 
tems with a 17(1) symmetry, the Ward-Takahashi iden¬ 
tity (WTI) provides an important constraint on the full 
vertex T^{K,K) which enters into the correlation func¬ 
tions. Given a mean field like self energy, it is possible 
to analytically solve the WTI, and obtain the full vertex 
function along with the full correlation function [20l |2I] . 

We define the generalized correlation function 


\G{K)r>^{K,K)G{Kh''{K,K)\ , (6) 


K 


where K = K + Q and j^{K,K) is a bare vertex. 
From this we have the density-density Xpp{Q) = P^^{Q) 
and current-current correlation functions Xjj{Q) = 
P^^Q),i,j & {1j 2,3}. The bare and full vertices sat¬ 
isfy respectively 

qpr{K,K) = G^\K)-Go\K), (7) 

qpT^^iK,K) = G-\K)-G-\K), ( 8 ) 


the transverse response P-^'^iQ) which contains no collec¬ 
tive modes: 


K 


^r(Q) =ETr| \G{Kh^iK,K)G{K) 


+ Fp^iK)Y{K,K)Fp^{K) 

- F,,{K)Y{K,K)F,c{K)] |>,( 10 ) 


where Fra{K) = AmGo{K)G{K) = F^{K) for m G 
{sc,pg}. Note that Fpg does not represent an anomalous 
Green’s function, but rather reflects a vertex correction 
to the correlation functions [5D] . There is a disagreement 
in the literature [nun] as to the importance of collective 
modes in the spin response below Tc- We agree with the 
results of Ref. [21], where collective modes were included. 

As shown in the supplementary material, when one 
integrates over the entire frequency range, a consequence 
of the WTI is that the /-sum rule is satisfied: 



with the latter a consequence of the WTI. We now spe¬ 
cialize to systems with the self energy as in Eq. ([^ . Using 
the WTI above Tc we have 

r'^(iF,iF) =7'^(i7,iF) + A2Go(W(iF,i7)Go(iF), (9) 

where ^^{K, K) = —K^ay is a time-reversed 

vertex. Inserting the full vertex into Eq. (§ then gives 
the correlation functions above Tc. 

One can incorporate superconducting (or equivalently 
superfluid) terms within this formalism building on 
Eq. (§ and, for example, address the superfluid den¬ 
sity |25| . as outlined in the supplement. One considers 


where X'k is the imaginary part of the density response 
function. This /-sum rule depends on the total particle 
number n and the bare mass m. Since A does not enter, 
the presence of spin-orbit coupling does not modify the 
weight of the /-sum rule. 

Spin response and f-sum rules.— In the spin channel, 
where there is no U{1) symmetry to justify the use of 
the WTI. Nevertheless, we are able to provide an a pos¬ 
teriori check on any proposed correlation function via a 
sum rule which we now derive. We define XsiS- (*<j-’) q) = 
J dr (Tr'S'qi (r) S-qj (0)) where T^ is the time order¬ 
ing operator and Sqi = (^i)ss’ Ck+qs' is the 

many-body spin density operator. Using the Heisenberg 
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equations of motion and the properties of Fourier trans¬ 
forms, the sum rule for the spin-spin correlation function 
X’g.g. can be shown to be 

/ V (^,q)) = ([[Ho,5q.] ,5_q,]), (12) 

where Uq = (k) Ck^' and Xg.s^ is the sin¬ 

gular part of XSiSj found by analytically continuing 
iuj ^ io + i5 and then taking the i5 —> 0 limit. 

Here we give the explicit result, for two example cases 
of interest and present further details in the supplemen¬ 
tary material: 

— {-uXs.s, (w; q)) =-V o;f^^nka■ (13) 

TT ^ ' m m 


where i € {x,z}, f^z = k±, and rika = 

TJ2iu, (^)i ■'^kh Gfj{K) a helicity Green’s functions 
to be defined in the next section. 

Correlation functions in the helicity basis.— In the ab¬ 
sence of a magnetic field, helicity is a good quantum 
number, and the correlation functions are most easily 
expressed in terms of the helicity Green’s functions [35] : 


Gf,{K) ^ . y 

lUJ -fi'ka 


Fh{K) = UkaVka 


ioj + Eka ’ 

1 


1 


iu) -I- Eka — Eka 


(14) 

(15) 


where Eka = \/fka + Cka = k'^/2m — ^ + aXk±/m is 
an eigenvalue of i?°(k), and Ff^{K) represents the pseu¬ 
dogap, or equivalently vertex contribution. Here o = ± 
denotes the helicity index and the coherence factors sat¬ 
isfy uL = + ^ka/Eka), mL + = 1 - 

It follows from the vertex function in Eq. ([^ that the 
explicit form for the /-sum rule | 20 j consistent density- 
density correlation function is 


Xpp(w,q) = i {1 + aa'cos {(j)k+q - M) 

K, 0 ',a' 

X 'Gfi{K)Gi{K) + FfiiK)F^'(K) . 


(16) 


The angle exp(i(/)k) = {kx+iky)/k±, so that exp(i(/)_k) = 

-exp(i())k)- 

The spin-spin correlation functions are constructed us¬ 
ing their form below Tc (deduced using the path inte¬ 
gral [ 21 ]) with appropriate sign changes in the pseudo¬ 
gap relative to the condensate gap. These sign changes, 
which appear in the T < Tf. Ward-Takahashi identity [37] 
are essential for satisfying sum rules. 

As can be shown, in the normal phase the following 
expression for the spin-spin correlation functions are fully 
compatible with the spin /-sum rules given in Eq. 


XSiS. (w, q) = ^ Y ± M) 

K,CK,(y' 

X G^(K)G“ (ix) + F^iK)E^\K) 


(17) 


where the -f, —signs are for Xs^^s^jXs^s^ respectively. 

Numerical Results.— We now look for qualitative new 
physics in the spin-spin response functions. We numeri¬ 
cally calculate the response function Xg g (q, w) at fixed 
q = (0.5, 0 , 0)kF as a function of w [3T], and for definite¬ 
ness consider T = 0.28Tf > Tc and unitary scattering, 
l/kpa = 0. We plot the results in Fig.l^ In order to illus¬ 
trate the physics, in Fig.j^a) and Fig.]^b), Rashba SOC 
or pseudogap effects were set to zero respectively, while 
Fig. [^c) shows their combined effects. The /-sum rules 
derived above are important for constraining numerical 
results of the spin-spin and density-density correlation 
functions. Gomparison between our numerical calcula¬ 
tions and the exact /-sum rules agreed to within a few 
percent. 

In Fig.j^a) we set A = 0. In this case, above Tc the spin 
and density correlations are equal {Xg^g^ = X'/p) and this 
function is plotted in the figure. Two low energy peaks 
are observable, as found in our earlier work [32] . The 
lower frequency peak reflects contributions from ther¬ 
mally excited fermions, while the higher frequency peak is 
associated with the contribution from broken pairs which 
appears at a threshold associated with the pseudogap. 

In Fig. j^b) we set Apg = 0 and plot Xg^g^ for a pure 
SOC system with A = 1.2kF. (We do not show X'/p since 
there is still no qualitative signature of A 7 ^ 0.) The 
response X'/ g shows two peaks, but one is at a consid¬ 
erably higher energy compared to Fig. [^a). The lower 
frequency peak reflects intra-helicity band contributions 
while the larger frequency peak is due to inter-helicity 
effects. 

Importantly, this figure shows how the physics of the 
Rashba ring band-structure can be directly probed by 
the spin-spin response function. To illustrate this, in 
the inset we plot the dispersion relation of two helicity 
bands. The horizontal line denotes the self-consistently 
determined chemical potential, chosen so that occupied 
fermions mostly reside in the Rashba ring. The onset of 
the inter-helicity band transition energy is given by the 
energy difference between two bands positioned on the 
inner circle of the ring, while the endpoint frequency for 
this peak is determined by the outer circle. These energy 
differences roughly match the width observed in the high 
frequency peak in the main plot. (The smearing of the 
width is because we have a non-zero momentum q and 
Tt^O.) 

Finally, in Fig. j^c) we plot X'g^g^ for the case where 
both pseudogap and Rashba SOG are present. Here 
we observe three distinct peaks. The first is associated 
with thermally excited fermions within the lowest helicity 
band, the second with the breaking of the preformed (pg) 
pairs and the third mainly with the inter-helicity tran¬ 
sitions discussed in the previous panel. We also observe 
some inter-play between pairing and the high frequency 
SOC peak, as this inter-helicity band peak is pushed to¬ 
ward slightly higher energies. 
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Conclusion.— A major finding of this paper is that 
spin-spin correlations provide a clear signature of the si¬ 
multaneous presence of Rashba modified band-structure 
and of a pairing gap. Signatures of both are a neces¬ 
sary (but clearly not sufficient) condition for ultimately 
obtaining a topological superfluid. This should comple¬ 
ment observations which are based on the single particle 
response functions in different experiments either in cold 
gases [TMTS] or in condensed matter. Our spin correla¬ 
tion functions are consistent with sum rules which we de¬ 
rive in this paper. These provide important constraints 
on the spin response which is complicated by the fact 
that spin conservation laws are unavailable for spin-orbit 
coupled systems. 
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Supplementary Material: Signatures of pairing and spin-orbit coupling in correlation 

functions of Fermi gases 

Chien-Te Wu, Brandon M. Anderson, Rufus Boyack and K. Levin 
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I. GREEN’S FUNCTIONS AND SELF ENERGY IN THE HELICITY BASIS 


In evaluating the correlation functions, it is convenient to have the Green’s functions written in the helicity basis. 
Here we present a brief discussion of the simple expressions used in the main text: 


G^K) ^ 

lUJ -£/ka 


Fh{K) = MkaWka 




lUJ - 
1 


A'ka 


1 


ioj + Aka iw — Aka 


( 1 ) 

( 2 ) 


where the Bogoliubov spectrum Aka = a/ with corresponding coherence factors which satisfy = |(1 + 
Cka/Aka)) Uka + ^ka = 1' Note that the helicity Green’s functions have the standard form, aside from the index a. 

We define the Nambu spinor = ^Ck^, Ck^, of annihilation and creation operators. The many-body 

Nambu Green’s functions are then defined as 


G{K) = - 




(3) 


(4) 


where P = I/T is the inverse temperature. The Green’s function G{K), time-reversed Green’s function G{K) = 
iay[G{—K)]'^iay, and the anomalous Green’s functions F{K) and F{K) are now in general 2x2 matrices. 

Using the equations of motion for the Nambu Green’s function in the pairing approximation we have 

^ ^ I A* G^\K), 

For the remainder of this supplement, as well as in the main text, we chose A to be real. The non-interacting 
Green’s function of the particle sector is Gq ^(A) = ioj — A°(k), and the hole Green’s function is defined through 
Gq^{K) = iaylG^^{—K)]'^iay = iuj + A°(k). This non-interacting Greens function depends on the single particle 
Hamiltonian through H^{h) = /2m — 71 + Akj_ • a /m where a = {a^, cfy, (Jz) is a vector of Pauli matrices in spin 

space. We consider a 2D Rashba SOG, so that a is coupled only to the in plane momentum kj^ = ky, 0). Taking 
the block matrix inverse gives the following 2x2 matrix Green’s functions: 

G(A) = {g^\K) + A2Go(A))”' , (5) 

F{K) = AGo(A)G(A). (6) 

In the pairing approximation, the structure of the block matrix inverse forces F{K) = A*/AF{K) = F{K). 

For pure SOG (i.e., no magnetic field), the system is time-reversal invariant and helicity is a good quantum number. 
Therefore, the hole Green’s function and the particle Green’s function can be diagonalized simultaneously with the 


same unitary transformation: Uk = exp 


with kj^ = kj^/fc_L and = {kx -f iky)/k_\_. This 


1^—i(/)k7!‘/4 ^kj_ X zj - a 

unitary transformation will then diagonalize the non-interacting Green’s function Go (A), and leads to the identity: 

I 


Go(A)=V . ^ Aka, 

^ ZW - ^ka 


(7) 


where Aka = Ak(l-l-acrz)C/j^/2 = (I-|-akj_ ■a/k±)/2 is a projector into the helicity state a = ±1. We can now express 
the G{K) and F{K) in terms of the projector Aka and the helicity Green’s functions as Eq. (I) and Eq. (2) through: 


G(A)=^G^(A)Aka, 

a 

A(A)=^A]^(A)Aka. 


( 8 ) 


(9) 
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II. THE WARD-TAKAHASHI IDENTITY 
A. Deriving the full vertex 

The Ward-Takahashi identity (WTI) in quantum field theory is a diagrammatic identity that imposes a symmetry 
between response functions. The particular symmetry we are interested in is the 17(1) abelian gauge symmetry that 
is present in the density/current channel. Since the pseudogap phase self energy is mean field, the WTI provides 
a powerful tool in deriving the full vertex that appears in correlation functions. This is because the WTI can be 
explicitly used to solve for the full vertex function, and thereby obtain the exact correlation function. 

We confine our attention here to the normal phase, as the density-density correlation function has collective mode 
effects which are complicated below T^. 

The WTI for the full vertex r^{K,K) is^ 

q^r>^iK,K)=G-\K)-G-\K). ( 10 ) 

Similarly the WTI for the bare vertex 7 ^( 77 , K) is K) = Gq^{K) — Gq^{K) and we also define K) = 

Gq^{K) — Gq^{K). Here K = K + Q. In the pseudogap phase the mean field like self energy is S(A) = A‘^Go{K). 
Using this self energy, the WTI becomes 

q^r^{K,K) = G^\k) - G^\K) - S(A) + S(A) 

= G^\k) - G^\K) - A2 (Go(A) - Go(A)) 

= q^r{k,K) - A^Goik) (g^\K) - Go 1(A)) Go(A). 

= q^nk, K) - A^Goik) [-q^Yik, K)) Go(A). 


It follows that the full vertex is then 

r^(A, K) = nk, K) + A^Goik)Y{k, A)Go(A). (11) 

It is important to note that, this expression is exact, and is not a perturbation expansion in terms of bare vertices. 


B. / and longitudinal sum rules 


We now show that, given the exact correlation function obtained by using the WTI, the /-sum rule is explicitly 
satisfied. This ensures our correlation functions are consistent with charge conservation. 

The correlation function is given by 


^ Tr [G(A)r^(A, A)G(A) 7 ^(A, K) 

K 

Inserting the full vertex from Eq. (11) gives the correlation function above Tc as: 

P^kQ) = E lGik)nk, K)G{K) + Fpg(A)7^(A, A)Fpg(A)]7^(A, A) 


K 


where Apg(A) = ApgGo(A)G(A), Upg(A) = ApgGo(A)G(A) = Apg(A). 

Applying the WTI to and using the cyclic property of the trace, then gives 

(g) - q • p- (g) = E G( A) [ 7 "(A, A + g) - (A - Q, A)] I . 

For the SOC system the bare vertex is 7 '^(A, A) = (l, -|- ^O’J.) ■ Using this yields 

HpO'^(g) - q • A'^(g) = ^(1 - So,.) E Tr {G(A)} = ^(1 - Jo,.)- 

m Z—/ rn 


( 12 ) 


(13) 


(14) 


K 


m 


( 15 ) 
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Here we have used n = Tr {G{K)}. In components, this equation becomes 

a;P°°(a;,q)-q-P"(u;,q) =0, (16) 

a;P°^-(cc,q)-q-P*^(a;,q) = -q. (17) 

771 

Setting a; = 0 in the second equation and then taking the dot product with q gives 

q-P*^( 0 ,q)-q = -^. (18) 

m 

Now use the identity ImP®°(a;, q) = —ImP°®(—w, —q) and solve for ImP°° in terms of ImP*-!. Applying the Kramers- 
Kronig relations and Eq. (18) then gives 


= -q-ReP*J (q,0) -q 

9 

nq^ 

m 

Defining Xpp{Q) = P^^{Q),Xjj{Q) = P^^{Q),i,j € {1j 2,3}, the /-sum rule, which holds for all q, is thus 



(19) 


( 20 ) 


Here the singular part of the correlation function, X"p(w,q), is equal to ImXpp(w,q), which is true only in the 
density/current channel. The next section on the spin-spin correlation functions elaborates further on this point. 
Similarly the longitudinal sum rule, which holds for all q in this continuum case, and above Tc, is thus 


/ 


du} ^ q-X 0 j(u;,q)-q ^ 


9 

nq 

m 


( 21 ) 


III. SPIN-SPIN CORRELATION FUNCTIONS AND SUM RULES 


A. General sum rules for spin density-spin density correlation functions 


Here we investigate both the spin-spin correlation functions and sum rules. In the presence of SOC, spin is no longer 
conserved. The WTI is therefore inapplicable to deduce the form of the spin-spin correlation functions. However, the 
Heisenberg equations of motion can be utilized to derive a sum rule for spin density-spin density correlation functions. 

We start with the definition for the spin-spin correlation function in the main text, XsiSj (*w, q) = / dr 
(PrS'qi (r) S'-qj (0)), for a many-body spin operator S'q^ = X^ss'k'^iLs (^i)ss' Ck-i-qs'- We will consider a general many- 
body Hamiltonian H — Ho +77/, where the contribution "Ho = X^ss'k '^ks-^ss' (^) is constructed from an arbitrary 
single particle Hamiltonian (k). We also assume that the interaction Hamiltonian commutes with the spin- 
operator, [Hi, S'qi] = 0. 

In order to constrain the spin-spin correlation function, we investigate the analogous sum rule for Xg.g_ {lo, q) that is 
derived in the previous section for the density/current response. To do this, we first analytically continue XsiSj (*^1 q) 
to real frequency by setting itu = uj + i6, then letting d —)■ 0, and finally invoking the Sokhotski-Plemelj theorem^. 
This is expressed mathematically as 


^SiSj (w, q) = lim^^o+S.Sj (7w, q) , 

iu)—co-\-i5 

= (w, q) + ix'ks, (w, q), 


= V 


dz q) 

TT Z — LO 


+ (‘k'; q)- 


( 22 ) 


Here V denotes the Cauchy principal value. Note that, for i ^ j, x!g.g_ (w, q) need not necessarily be real^. Returning 
to the time domain, the singular response function is X'/.g. {t, q) = 5 ([•S'qi (t ), S-qj ( 0 )]). 
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We now can derive the sum rule. Using the properties of Fourier transforms, we have 

J ^ ) 

= ([[H,^q,],5_q,]), (23) 

where we have used the Heisenberg equations of motion idtSqi (t) = — [H, Sqi{t)] and the commutators on the last 
line are taken at equal time. The problem thus reduces to calculating the equal time commutator {[[H, •^qi] , >5'-qj]) — 
["Ho, S'qi] . 

The commutator in Eq. (23) is tedious but straightforward to calculate. After much algebra we arrive at the /-sum 
rule, or analogue thereof, for spin density-spin density response: 

/ ^ (w,q)) = (k) -i?°(k + q)}CTj -CTj {cr*,i?° (k-q) - i^°(k)} 

+ [Ef"(k)+ffO(k + q),a,]a,-a, [i/°(k-q) + /f«(k),a,]') (c^cus') ■ (24) 

This expression can be further simplified as follows. The Hamiltonian can be expressed as (k) = (k) where 

(k) = ^Tr [ui_iH^ (k)] and = (l,cr). After using the commutation and anti-commutation rules for Pauli matrices, 
further algebra reduces the sum rule to the form 

/ V (w, q)) = [iH° (k + q) + (k-q) - 2H° (k)) 5 k] S,, 

- 2 ((h (k + q) -b h (k - q)) • Sk% -b {hi (k + q) -b hi (k - q)) Skj) 

- i [{hk (k -b q) - hk (k - q)) rik - (/lo (k -b q) - /lo (k - q)) s^k] e*jfc|, (25) 

where we have defined g\^=T J2iuj G{iuj, k), nk = Tr [^k] and Ski = Tr [ui^k]- The first line, appearing only for i = j, 
can be seen to be equal to the density /-sum rule contribution. The second line can give diagonal contributions, 
resulting in a deviation of the spin density response from the density response in the presence of SOC. The third line 
vanishes for i = j. It is of note that this line is purely imaginary, and can give an imaginary /-sum contribution. This 
is only possible for correlation functions of two distinct operators, as was noted in the literature^. This implies that 
X'gig. 7 ^ ImXs-s^ for i 7 ^ j. Finally, the first and third lines of Eq. (25) are symmetric under interchange of f -f-)- / and 
q e7 —q, while the second line is not. This can be explained by observing that Eq. (23) is not fully symmetric under 
the the exchange of i e 7 / and q e 7 —q. 


B. Explicit spin-spin correlation functions and sum rules for Rashba SOC 

We now present simplified expressions for the correlation functions and sum rules for the case of Rashba SOC. Due 
to symmetry in the x-y plane, there are only four unique correlation functions of interest. The remaining spin density 
correlation functions can be found through exchanging x ^ y where relevant. The explicit spin-spin correlation 
functions are 


^S,,S,,(w,q) 
(w,q) 
Xs,S, (w,q) 


E (l + aa'cos(0k+q + <(>k)) \G%{K)Gi{K) + F^{K)FS' {K) 


K,oc,oc' 


I ^ (aa'sin(</)k+q) - sin(0k)) [g^(K)G“ (K) + F^{K)F^ {K) 


K^Ol^Oc' 


1 ^ (1 _ cos (0k+q - <(>k)) [G^(i^)G“ {K) + F^{K)FS {K) 


Ka 


Xs,s„(w,q) = 0. 


(26) 

(27) 

(28) 
(29) 


Note that the response Xs^Sy vanishes identically, as it must due to in-plane rotational symmetry . 
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For each of the above correlation functions, we have explicitly integrated the left hand side of the sum rule, Eq. 
(25), and found: 



where n^a = Tr [Pkaffk] = TJ2iuj These results are consistent with the form using the right hand side of Eq. 

(25). Note further that the sum rule for Xs^S;, is purely imaginary. This is not surprising, as mentioned above. 


IV. SUPERFLUID DENSITY 


The formalism underlying this paper is capable of addressing the superfluid phase. However, many response 
functions in this phase require the inclusion of collective mode effects. The superfluid density, which can be obtained 
using a transverse response, does not require collective modes. The superfluid density is defined by 




(33) 


where XjjiQ) is the transverse part of the current-current correlation function. 

This contains two contributions from: (z) quasi-particles which are reflected in the current-current correlation 
function and (ii) the diamagnetic current which above Tc must precisely cancel the quasi-particle term. Below Tc 
the current-current correlation function acquires an addition contribution which depends on the anomalous Green’s 
function Fsc- Note that the diamagnetic current is dependent on the total pairing gap. The current-current correlation 
function depends on the function Fpg as defined above, and arises from vertex corrections. 

Using Eq. (10) from the main paper, the current-current correlation function and diamagnetic response tensors are 


tjj{Q)=J2^T 

K 

m / _ / 


dia 


K 


J(k, q)G(iF)J(k, q)G(iF) + J(k, q)Fpg(i^)J(k, q)Fpg(iF) - J(k, q)Ese(iF)J(k, q)Ese(iF) 
J(k, 0)G(iF)J(k, O)G(iF) + J(k, 0)Epg(iF)J(k, O)Epg(iF) -b J(k, 0)F,c(iF)J(k, Q)F,,{K) 


(34) 

(35) 


where J(k,q) = is the current operator. Notice that the relative sign between the Fgc and Fpg 

contributions is different in XjjiQ) in the diamagnetic current. 

The superfluid density tensor is thus 

/<-> \ 


— =2^Tr J,(k,0)E,c(iF)Jj(k,0)F,e(iF) 


K 


(36) 


Due to symmetry, the off-diagonal (z 7^ j) elements vanish. Furthermore, the SOC leads to different in-plane {i = j = 
X, y) and out-of-plane (i = j = z) response. Despite a significantly different form, we can show this expression for the 
superfluid density is consistent with that derived by adding a phase twist to the order parameter^^®. 
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